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Abstract 

A semi-microscopic model to study the neutron and proton induced backbending phenomena in 
some deformed even-even nuclei from the rare earth region, is proposed. The space of particle-core 
states is defined by the angular momentum projection of a quadrupole deformed product state. The 
backbending phenomena are described by mixing four rotational bands, defined by a set of angular 
momentum projected states, and a model Hamiltonian describing a set of paired particles moving in 
a deformed mean field and interacting with a phenomenological deformed core. The ground band 
corresponds to the configuration where all particles are paired while the other rotational bands 
are built on one neutron or /and one proton broken pair. Four rare earth even-even nuclei which 
present the second anomaly in the observed moments of inertia are successfully treated within the 
proposed model. 

PACS numbers: 21.10.Re, 21.60.Ev, 21.10.Hw, 27.70. +q 
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I. INTRODUCTION 



The irregular behavior of the moment of inertia in the yrast band at intermediate and 
high spin states, known as backbending, has always attracted considerable experimental and 
theoretical attention. Since its first experimental observation many endeavor attempts 
were performed in order to explain the phenomenon. It is commonly accepted that it is 
caused by the intersection of two rotational bands. This interpretation was proposed by 
Stephens and Simon |2| based on the rotational alignment of the individual single-particle 
angular momenta of a broken pair along the rotation axis. The pair breaking is caused by the 
Coriolis force which violates the time-reversal symmetry. The first theoretical interpretation 
based on the Coriolis anti-pairing effect was due to Mottelson and Valatin {^J where the 
backbending phenomenon was put on the account of a drastic change in the pairing field. 

Although the band hybridization method was all along known and applied to this par- 
ticular problem within some phenomenological approaches 4|-|6|], the nature of the involved 
rotational bands was not yet well established. Only after the rotational alignment hypothe- 
sis was confirmed, it became clear that the first backbending is due to the intersection of the 
ground band (g) and a two quasiparticle (2qp) band built upon a broken pair from a high 
angular momentum orbital. The second band is often referred to as the S'(tockholm) band. 
Thus, the anomalous increase of the moment of inertia is interpreted as the reduction of the 
energy cost to achieve a certain total angular momentum by aligning the angular momenta 
carried by the constituents of a broken pair. Stephens and Simon noticed that in the rare 
earth region the first broken pair is from the neutron intruder orbital 62 13 / 2 . Actually this 
picture was later confirmed by many theoretical calculations, mostly based on the cranking 
Tartree-Fock-Bogoliubov (CHFB) p], |8[ calculations and the core plus quasiparticle models 
, [loj]. The backbending is a relatively widespread phenomenon within the rare earth 
region, but only very few nuclei exhibit a second anomaly in the moment of inertia. It was 
for the first time measured for 158 Er [111], and the early interpretation was based on the 
alig nment of the individual angular momenta resulting from breaking a 5/iii/2 proton pair 



121 ] . Other nuclei which exhibit a second moment of inertia anomaly are located around 



the N = 90 rare earth isotopes. The proton nature of the second broken pair was at a first 



glance queried in Ref.[13|, since in the same energy region of the spectrum, the alignment 
of a 5hg/ 2 neutron broken pair might also play an important role. However, the proton 
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nature of the second backbending was later confirmed by several more detailed theoretical 
studies j3, Jj| based on blocking arguments offered by the experimental investigations of 
the odd-proton and odd- neutron neighboring nuclei of the N m 90 isotopes jl^, As a 
result, the second backbending is regarded as being caused by a successive breaking of a 
neutron and a proton pair, where the neutron broken pair is the one which causes the first 
backbending 17J. As a matter of fact, the suspected neutron pair 5/19/2 which may break at 
a time with 5hu/2 proton pair is causing, indeed, a third anomaly in the moment of inertia 
of some isotopes of Yb [18]. Indeed, for this nucleus a weak up-bending is noticed at spins 
beyond J = 36. 

Of course band hybridization is a conventional name used in the early publications quoted 
above. Hereafter we replace it by band mixing which might be used in a brighter context. 

The most extensive calculations on the double backbending were performed in the frame- 
work of the CHFB approach, which provided one of the most reliable qualitative description 
of the phenomenon over a large number of nuclei. One of the most important features of the 
CHFB approach is that it embraces all the mechanisms known to cause the backbending, 
that is the particle alignment, the pairing phase transition and the sudden change of defor- 
mation. However, the CHFB description is a semiclassical one, which encounters difficulties 
in describing the states near the band crossing. An important improvement is obtained 
by the angular- momentum-projected Tamm-Dancoff approximation which was successfully 
applied for the dysprosium isotopes [li], Q. Therefore, in order to achieve a quantitative 
description of the multiple backbending, a full quantal formalism is necessary. Such models 
were proposed based on mainly two directions: genuine shell model formalisms [n| can trace 
better the influence of the single-particle degrees of freedom on the pair breaking process 



while the particle-core models 
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21] put emphasis on the rotational alignment description. 



22 



23J can be also included in the 



The calculations based on the interacting boson model 
first category. For a quantitative description of the energy spectra with double backbending 
one advocates for the second solution. The advantage of the particle-core approach con- 
sists in the fact that it treats the single-particle and collective degrees of freedom on equal 
footing. It is worth mentioning that a qualitative explanation of the first backbending in 
some isotopes of Pt, W and Os, was obtained in Ref.|24j by using the general collective 



model 
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261 ] where, of course, the particle degrees of freedom are missing. Therein the 



backbending is determined by the angular momentum dependence of the moment of inertia, 
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induced by the specific ways the structure coefficients are fixed. 

In a previous publication [27], we proposed a semi- microscopic model for the description of 
the backbending phenomenon within the band mixing picture. The rotational bands implied 
in the mixing procedure were defined by angular momentum projection from quadrupole 
deformed product states and a model Hamiltonian describing a set of intruder neutrons 
interacting among themselves through pairing forces and coupled to a phenomenological 
deformed core. By projecting the angular momentum one avoids the difficulties showing 
up when one treats observables which are sensible to the angular momentum fluctuations. 
Indeed, working with states of good angular momentum is more advantageous than applying 
cranking methods which encounter enormous angular momentum fluctuations in the band 
crossing region. The distinctive feature of our model is that, although we use a spherical 
projected particle-core basis, the core and the single-particle trajectories are deformed. The 
mixing of the rotational bands was achieved by diagonalizing the model Hamiltonian in an 
orthogonal basis constructed from the projected states of g and S"-bands. The model was 
meant to reproduce only the first backbending, which was done quite well for six even-even 
nuclei from the rare earth region. Besides the reproduction of the backbending plots, the 
formalism {27] also provided some useful information regarding the rotational alignment of 
the particles moving in an intruder orbital. 



In the present paper we extend the formalism from Ref. |27j to the second backbending 
induced by a proton broken pair. This is done by performing the mixing of four rotational 
bands. The first two are obviously the g-band and the S'-band with a neutron broken pair, 
whereas the other two are associated to a proton broken pair and to two, one of neutron and 
one of proton type, broken pairs, respectively. The projected states which define the four 
bands have specific single-particle factors describing each case mentioned above. The protons 
and neutrons are treated through BCS model states associated only to 6z/i 13 / 2 and 57r/in/2 
orbitals. The intruder particles are coupled to a phenomenological core which is deformed 
and described by means of the coherent state model (CSM) [28|]. The projected states are 
deformed and therefore not orthogonal but can be used to construct an orthogonal basis. 
The lowest eigenvalues of the model Hamiltonian in this orthogonal basis define the yrast 
band. The main purpose of the present work is to reproduce the experimental yrast spectrum 
and its backbending behavior for some even-even rare earth nuclei which are known to be 
double backbenders, as well as to provide a through out analysis of the rotational alignment 
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process and the possible consequences for the E2 transition properties along the yrast band. 

The description of the method and results are presented according to the following plan. 
The model Hamiltonian as well as the projected particle-core product basis used for the 
description of the double backbending phenomenon is presented in the next section, Section 
II. The E2 transition probabilities are considered in Sec. Ill and the emerging numerical 
calculations are given in Sec. IV. Final conclusions are drawn in Sec. V. 



II. THE MODEL HAMILTONIAN AND A PARTICLE-CORE PRODUCT BASIS 

In this study we present a new and simple semi-phenomenological model to explain the 
first two backbendings seen in some rare earth nuclei. The spectra exhibiting a double back- 
bending will be described by a particle-core Hamiltonian whose eigenvalues are calculated 
within a particle-core space. The nucleons are moving in a deformed mean field and the 
alike ones interact among themselves by pairing force. The core is deformed and described 
by a phenomenological quadrupole coherent state [281 ]: 

Vc = e d ^o-M| } 6) (2.1) 

where fejL with —2 < \i < 2 denotes the quadrupole boson operator, while d is a real 
parameter which simulates the nuclear deformation. The two subsystems interact with each 
other by a qQ and a spin-spin, Jf ■ J c , interaction. The associated Hamiltonian is: 

H = H c + H sp + Hpair + H pc . (2.2) 

The core term H c is a quadratic polynomial of the quadrupole boson number operator, 

H c = u b N + uj\N 2 . (2.3) 

As for the single-particle Hamiltonian H sp , this is a sum of two terms corresponding to 
neutrons and protons, each of them describing a set of particles in an intruder spherical shell 
model orbital \nlj): 

Hsp = { e njkji — \) ^mlijim^nikjimi- (2-4) 

i=i/,Tr mi=all 

Here - m and c n ij m are the creation and annihilation operators for a particle in a spherical 
shell model state \nljm) with the energy e n ij, while A is the Fermi level energy for the system 
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of paired particles. Alike nucleons interact through a pairing force: 



Hpair — ^ ^ PjiPjn (2-5) 

i=V,1T 

where Pj and Pj denote the creation and annihilation operators of a Cooper pair in the 
intruder orbital j. 

The particle-core interaction consists of two terms, the quadrupole-quadrupole (qQ) and 
the spin-spin interaction: 

Hpc = HqQ + Hj f j c , 

H qQ = -A c Yl (^ / ^ m *l r2y 2Ml^^X)^^m l c ^^K (-)" 6 2-„ + K > 
H JfJc = CJ f ■ J c . (2.6) 
Here the total angular momentum carried by protons and neutrons is denoted by: 

Jf = J p + Jn- (2.7) 

The interaction strength Ac is taken to be the same for neutrons and protons. The pa- 
rameters Ac and C are free in the present work and therefore are to be fixed by a fitting 
procedure. 

The mean field is defined by averaging H(= H sp + H q q) with the coherent state ( 12.11) . 
which results in obtaining a single-particle Hamiltonian which is similar to the deformed 



Nilsson Hamiltonian 



34|. In the first order of perturbation, the energies of the deformed 



mean field are given by: 



£nijm = e„ij - AdX c (2n + 3)C{ l{C J m 2 J w with 

2 u 2 



x ° = sL f* Ac ' (2 - 8) 

where n is the principal quantum number of the intruder orbital, while M and Uq are the 
nucleon mass and the harmonic oscillator frequency. e„.u denotes the spherical shell model 



energies corresponding to the parameters given in Ref. 



351 ] i.e., 



ftwo = 41^T 1/3 , C = -2hw K, D = -hu KfM. (2.9) 

where the parameters have the values (0.0637,0.42) for neutrons and (0.0637,0.6) for 

protons. The second term in the right hand side of Eq. (l2.8j) is obtained by averaging the 



non-spherical part of the mean field with the spherical shell model state \nljm). The true 
eigenvalues of the mean-field would be obtained by diagonalization, when the off-diagonal 
matrix elements of the deformed term are taken into account. 

Of course one could argue that the single-particle energies with linear dependence on 
the deformation, look unrealistic. One undesired feature is that the state with j = 1/2 is 
not affected by deformation. Actually, we were aware of this drawback and corrected for 
it Q, 0- 

Briefly, a quadratic term in d could be obtained for example by adding the 
second order perturbative correction or by adding the monopole-monopole interaction to 
the particle-core Hamiltonian and then applying the first order perturbation theory to the 
unperturbed spherical term, as we actually did in the quoted references. Diagonalizing the 
mean-field Hamiltonian in the " asymptotic" basis one obtains the Nilsson energies and wave 
functions. 

We opted for the linear dependence on d for energies and the deformed basis \nljm) 
because of advantage of having the angular momentum as good quantum number and that 
happens despite the fact that the states are deformed. The mentioned problem of j = 1/2 
does not matter at all here since the intruders have high angular momenta. Moreover, for 
small deformation the single-particle energies approximate reasonably well the Nilsson ones. 
The pragmatic feature which is worth to be mentioned refers to the fact the the Fermi 
level in our model corresponds to the sub-state m which is equal to the Q associated to the 
Fermi level from the Nilsson scheme. This feature gave us the certainty that the essential 
ingredient for approaching the backbending behavior, is included. 

Pairing correlations with such a deformed basis but in a different context has been also 
used in Ref.J^sj]. Since only the relative energies to the Fermi level are involved in the BCS 
equations, the orbital energy e n ij is taken to be zero. Moreover, due to the fact that the 
quantum numbers n and / do not change within a multiplet we simplify the notation and 
denote the resulting energies by Sj m . From here it is obvious that two states related by a 
time-reversal transformation have the same energy, and therefore one can restrict the single- 
particle space to the states \jm) with m > 0, keeping in mind that each such state is occupied 
by a pair of nucleons. The sum of the mean field term and the pairing interaction for alike 
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nucleons is brought to a diagonal form through the Bogoliubov-Valatin (BV) transformation: 

«i* = u jk c) k - v jk (-y- k Cj _ k , 

a jk = U jk c jk - V jk (-y- k c]_ k . (2.10) 

The output of the BCS calculation consists of the occupation probabilities of the m- 
substates, the gap parameter A, as well as the Fermi energy A. Consequently the average 
number of nucleons in the j-multiplet, 2(N^ ir ), with r = v>, n is readily obtained: 

(N T p L) = J2 V im- (2-11) 

m>0 

For the chosen nuclei the Fermi levels for neutrons and protons, lie close to a sub-state of the 
intruders 6*13/2 and 5h u / 2 respectively. If the particle-core basis was a deformed one, then 
the lowest state \2qp)\ip c ) would correspond to a sub-state of the two intruders, respectively. 
The mentioned substates have m— 1/2 for neutrons and m = 7/2 for protons. Since the core 
state does not contribute to the total K quantum number, the projection of the total angular 
momentum on the symmetry axis, we say that the intrinsic states leading to the yrast band 
have a K = 1/2 for the neutrons and K = 7/2 for protons. Also in the Nilsson model, 
the last filled neutron state has = 1/2 while the last proton occupies the state Q = 7/2. 
The choice of the K = 1/2 sub-state as the Fermi level of the neutron system was made in 



Ref. [27| to describe the first backbending. As for considering the K = 7/2 Fermi level for 
the proton system, breaking the corresponding pair and aligning the resulting quasiparticle 
angular momenta to that of the core as prerequisite conditions of the second backbending, 
these features are in full agreement with the microscopic formalism used in literature. In 
this respect in Ref.Q] the alignment of a = 7/2 broken pair is used to explain the second 
backbending in 158 Er and 160 Yb. The last nucleus mentioned is also treated by Cwiok and 



collaborators in Ref. [29] while 158 Er by Riley 



15| 



A great simplification is obtained if the single-particle space is restricted to the intruder 
multiplets where a number of nucleons equal to 2{N^ 3 air ) is distributed. Solving the BCS 
equations in the restricted space, the quasiparticle energies depend on m but are still in- 
variant at changing m to —m. However, in a pure microscopic formalism where the Coriolis 
interaction is included in the mean field, the time reversed quasiparticle states are no longer 
degenerate and consequently the broken pair is a K = 1 state. Here the term Jf ■ J c , which 



simulates the Coriolis interaction in the sense specified in Ref. 



271 ] is only subsequently 



used, when the whole Hamiltonian is diagonalized and thereby the broken pairs with K = 1 
are used. An important technical simplification is achieved if these pairs are obtained by 
applying the angular momentum raising operator on the K — pairs. 

If the quasiparticles were not deformed and moreover the dangerous graphs were elimi- 
nated at the level of BCS calculations, one would expect that the interaction between states 
with different number of quasiparticles is vanishing. Under these circumstances, truncating 
the particle-core space to the states with Oqp, 2qp and 4qp is a reasonable approximation. 
Since the rotation process involved in the angular momentum projection operation changes 
the K quantum number, and moreover particles and holes are mixed by the BV transfor- 
mation, the overlap of states with different number of particles is however nonvanishing. 
Despite this feature we keep the restriction of the quasiparticle space as specified above. 
The reason is that the mixing weight of components with more than 4 quasiparticles would 
be at least of sixth order in the U and V coefficients and consequently small. 

Thus, the restricted space of angular momentum projected states to be used for treating 
the model Hamiltonian, H qp , written in the quasiparticle representation is: 

{*Sf. *ji;l(j>), ■ (2.12) 

The set members are defined by: 

= Mf ] PUnBCS) d \pBCS) d ^ (2.13) 
^ 2 i ; i0» = A4?0»P^ [j + a) nU a) n _ v \nBCS) d ] \pBCS)^ c , (2.14) 
= Nfl{ Jp Tt)P( n \nBCS) d \j + a] v AApBCS) d ] ^ (2.15) 



^JM:20>;j>) = ^ ( SUn^3p^)P J M2 J +®} n „ a ] n -v\ nBCS ) d J+a}„*<Xj p -n \pBCS) 



(2.16) 



where the reciprocal norms can be analytically expressed. Also, at /o^ stand for the 
creation/annihilation quasiparticle operators. The Hill- Wheeler projection operator [30j has 
the form: 

Pmk = I D{^ K R(Q)dQ. (2.17) 



8tt 2 

The angular momentum projection from the many body fermion states is achieved by using 



the procedure of Ref. 



3 if ] - The Pauli principle restrains the maximal angular momentum of 
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a given configuration |32| to 

Jr x = - 2iV p T ^ + 1), (2.18) 

where Nl^ ir pairs of r particles, occupy the states of angular momentum j T . 

The set of projected states mentioned above, is not orthogonal. We orthogonalized first 
the angular momentum projected basis and then diagonalized the model Hamiltonian written 
in the quasiparticle representation. Note that the bands mixing is achieved by two processes, 
the orthogonalization procedure of the initial basis and then by diagonalizing the model 
Hamiltonian H (I2.2p . The lowest eigenvalues of the total Hamiltonian H in the orthogonal 
basis defines the yrast band. 

The energy spectrum of the rotational bands is approximated by the average of the total 
Hamiltonian with each projected state from the set 02.121) . 

n 

The mixing of these bands is achieved following the procedure of Ref . [27| extended to 
the case of four interacting bands. Here we briefly present the main ingredients of this 
procedure. 

Indeed, denoting by a J m the eigenvalues and by V? m the eigenvectors of the overlap matrix 
corresponding to J ^ 0, it can be checked that the set of functions 

Y,^mVL m = l,2,3,4, (2.19) 



V "m i=l 
is orthogonal. 

Writing the total wave function as an expansion in the newly obtained orthogonal basis: 

^ = X>Xf, (2-20) 

m=l 

the eigenvalue equation associated to the model Hamiltonian acquires the following matrix 
form: 

4 

E ^' X -' = E ™ XJ ™- ( 2 - 21 ) 
m'=l 

The Hamiltonian matrix is defined as 



i 



Kl> = ^F==j= ^ V^Z\H\^)V^, m , (2.22) 

Solving the homogeneous system of linear equations (I2.2ip for a given J ^ and then 
changing J, one obtains a four J-sets of energies. Collecting the lowest energy from each 
J-set of solutions, one obtains the so called yrast band. 
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III. E2 TRANSITION PROBABILITIES 



The reduced quadrupole transition probabilities are calculated by truncating the transi- 
tion operator to the boson part, i.e. we suppose that the collective transition is due to the 
core component of the wave function. The microscopic structure of the yrast states have 
however an indirect contribution. The boson structure of the transition operator is assumed 
to be of the form: 

= q[a 2fl + q' 2 (aa) 2fl , (3-1) 

where a 2fJi denotes the quadrupole coordinate which is depending linearly on the boson 
operators 

a 2lM = -L(bl + (-rb 2 -,). (3.2) 
In terms of quadrupole bosons the transition operator has the expression: 

Q2, = qi (4„ + (-Th,-,) + q2 + 2(4^ + {b 2 b 2 ) 2 ,) , (3.3) 

where: 

qi = -j=<&, qi = \q*2- (3-4) 

The reduced probability for the quadrupole transition in the yrast band, using the Rose's 



convention 



33|, can be written as 



2 



B(E2,J+ J'+)= (® J Tot \\Q2W Tot ) , (3.5) 

where the functions involved are the states (I2.20p obtained by diagonalizing the matrix ■ 
If the final state is + , then instead of $^' ot with J' = we use . The transition matrix 
elements involve two parameters q\ and q 2 , which are to be fixed by a fitting procedure. The 
reduced matrix elements of the transition operator have been analytically expressed in Refs. 



3, 



50|. 



IV. NUMERICAL APPLICATION AND DISCUSSIONS 

There are very few nuclei in the rare earth region which present a second anomaly in 
their moment of inertia evolution along the yrast band. The most studied nuclei are 156 Er, 
158 Er, 160 Yb and 162 Hf since for them a great deal of experimental data are available. These 
nuclei will be treated within the formalism described in the previous sections. 
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A. Parameters 



The model involves seven parameters. Six of them, namely the neutron and proton 
pairing constants G n and G p , the strengths of the qQ and spin-spin interactions, X c and C, 
and the strengths Uq and w\, of the two boson terms, are the structure coefficients defining 
the model Hamiltonian. The remaining parameter d defines the coherent state ip c and plays 
the role of the deformation parameter. The fitted values of these parameters are given in 
Table I. In what follows we shall explain how these parameters were fixed. 



TABLE I: The fitted parameters for the four nuclei are listed. The nuclear quadrupole deformation 



02, taken from Ref. 



381 ] , is presented for comparison with the deformation parameter d. 



Nucleus 


d 


X c [keV] 


G n [MeV] 


G p [MeV] 


w§ [MeV] 


uj\ [keV] 


C [keV] 


d ■ X c [keV] 


P2 


156 Er 


1.9498 


84.0455 


0.2146 


0.2626 


1.1420 


0.255 


3.042 


163.87 


0.177 


158 Er 


2.4910 


68.6731 


0.1803 


0.2593 


1.1525 


-1.426 


5.866 


171.06 


0.203 


160 Yb 


2.2870 


74.9940 


0.1892 


0.2619 


1.2684 


-0.514 


2.270 


171.51 


0.195 


162 Hf 


2.1490 


78.2942 


0.2000 


0.2583 


1.3104 


8.674 


-1.991 


168.25 


0.184 



In the first step, the BCS equations were separately solved for protons and neutrons. 
The pairing constants and the single-particle energies represent the input data for the BCS 
equations. The single-particle energies are defined by Eq.f l2.8j) and depend linearly on the 
deformation parameter, as can be seen from FigJT] and FigfSJ From these one can see that 
the product dXc plays the role of the deformed mean field strength, like the quadrupole 
nuclear deformation fl 2 in the Nilsson model 34|. Given the fact that here we deal only with 
neutrons from the 6^13/2 intruder orbital and protons from 5hu/2 intruder orbital, which are 
responsible for the first and the second band crossing respectively, the BCS equations are 
solved only for a subset of the entire neutron and proton single-particle space which contains 
the states that might interact with the mentioned intruder states. Since the single-particle 
energies yielded by the deformed mean-field are m-dependent quantities, the substates of 
the intruders will be specified by adding a lower index m to the standard notation specific 
to the spherical single-particle states. Thus, for neutrons, the subset comprises all states 
of the n = 5 shell, excepting the substates with \m\ < 11/2 coming from 5/tn/2 orbital, 
together with the intruder states 6213/2,™ and the state 5/in/ 2 ,n/2 coming from below, which 
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is an intruder for the n = 4 shell. Similarly, the proton subset comprises all states of the 
n = 4 shell, the intruder state 4gg/ 2i 9/2 for the n = 3 shell coming from below and of course 
all intruder states 5/in/2,m- in total, one has to solve the BCS equations in a space of 23 
neutron states and 17 proton states where each single-particle state can accommodate two 
nucleons. The nuclei 158 Er, 160 Yb and 162 Hf are N = 90 isotones, such that we distributed in 
the neutron subspace 10 particles for each, and 20, 22 and 24 particles in the proton subspace 
respectively. As for 156 Er, this has 8 neutrons and 20 protons distributed in the corresponding 
subspaces. Judging from the observed degree of the shell filling, for all considered nuclei the 
last occupied proton intruder state has the projection 7/2, while the neutron intruder 
state i\3/2 which is closest to the neutron Fermi level has the projection 1/2. Thus, the m 
substates which correspond to the broken neutron and proton pairs (m^m^) = (1/2,7/2) 
is the same for all four nuclei, even though they have different neutron and proton numbers. 

The pairing interaction constants G n and G p and the qQ interaction strength are fixed so 
that the observed sequence of the single-particle levels and the last occupied state for a given 
deformation d of the core are reproduced. Later on, a fine tuning is performed in order to 
improve the position of the band crossing points. Solving the BCS equations one obtains the 
quasiparticle energies, the gap parameter A, the Fermi level energy A and the occupation 
probability parameters U and V. The projected neutron and proton single-particle states 
(I2.13p -f l2.16p describe only the nucleons from the intruder orbitals 6^13/2 and 57r/ixi/2- Thus, 
in further calculation one would need only the BCS parameters concerning the seven neutron 
states ix3/2 and the six proton states hu/2- Using the occupation probabilities of the intruder 
states, one calculates the average number of pairs in the considered intruder orbitals, {N^ 3 air ). 

It is needless to say that the BCS calculations performed only for the single-particle states 
of the considered intruder orbitals with a number (N^ ir ) of occupying pairs would provide 
results equivalent to those obtained for the larger single-particle subspaces chosen above. 
Even though the equation ( 12. 18ft is designed for an even and integer number of pairs, it can 
be used to determine an approximate higher limit of the angular momentum realized in a 
virtual configuration of (N^ ir ) pairs. The value obtained in this manner is then rounded 
to the closest even integer, defining in this way the upper limits of the summations over 
neutron and proton angular momenta J n and J p involved in the definition of the projected 
single-particle states. All this information and the BCS results are given in Table II. With 
all these data, the projected states (12. 131) - (12. 161) are fully determined. 
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Concluding, the BCS calculation in the extended single-particle space is used to calculate 
the average number of pairs in the intruder orbitals. Once these are determined we solved 
the BCS equation for each intruder keeping the obtained restriction for the number of the 
r-particles. Also using the average number of the r-pairs we calculate the maximal values 
of the angular momentum carried by the given system of fermions. In this way the space of 
the four particle-core projected states is readily defined. Thus, we stress again the fact that 
our method is based on a single j calculation and not on a many j. We used the many j 
calculation just to remove the ambiguity in determining the number of the r-nucleons which 
should be distributed among the intruder substates. 

TABLE II: The neutron and proton Fermi level energies, gap parameters and the quasiparticle 
energies are given for each treated nucleus. The average number of pairs determined with (|2.11|) 



and the corresponding exact and approximated maximal angular momenta (|2. 18|) . obtained by 
replacing the number of pairs Nl^ ir by the average number (Nl^ ir ), are also given. 
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i56 Er 


48.350 


1.39475 


1.44662 


0.95 


11.51 


12 


44.271 


1.46021 


1.46081 


3.33 


17.79 


18 


158 Er 


48.496 


1.13589 


1.13905 


1.26 


14.50 


14 


44.083 


1.39234 


1.39288 


3.34 


17.77 


18 


160 Yb 


48.268 


1.22848 


1.23084 


1.29 


14.74 


14 


44.395 


1.35984 


1.46115 


3.71 


16.98 


16 


162 Hf 


48.056 


1.34656 


1.34681 


1.30 


14.78 


14 


44.678 


1.27407 


1.61959 


4.05 


15.81 
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Some remarks concerning the BCS results are worth to be made. The observed single- 
neutron level structure of all four nuclei for the tabulated values of the quadrupole defor- 
mation /3 2 shows that none of the neutron intruder states i 13 / 2 are occupied. In the present 
model, the single-particle energies ( 12. 8 j) depend linearly on the deformation parameter d 
contrary to the Nilsson case. Because of this feature, one finds that for the N = 90 isotones 
the Fermi level provided by the BCS equations is right above the first intruder state 61^13/2, 
as indicated in figures [I] and EJ Exception is for the 156 Er isotope which has fewer neutrons 
and cannot fill any intruder 61^13/2 sub-state, but due to the large value of the pairing 
strength G n the occupation probability is considerably extended above the Fermi level A n 
and thus placing an average number of two nucleons in the intruder orbital 213/2 (see Table 
II). Also, according to the shell filling, the proton Fermi level of the Er isotopes must be 
placed under the intruder state 57r/in/ 2 ,7/2> but as can be seen in FigJTJthe obtained Fermi 
level A p is placed right above this state. This is caused by the fact that the intruder state 
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FIG. 1: Neutron (top) and proton (bottom) single-particle energy levels given in units of 
huj (= 41 A' 1 / 3 MeV) and calculated with Eq.flZSD for 156 Er (left) with X c = 84.0455 keV and 
158 Er (right) with Xc = 68.6731 keV. The vertical lines indicate the single-particle configurations 
corresponding to the fitted deformation parameter d. The Fermi energy level resulting form the 
BCS calculations is also pointed out. 

/ill/2,7/2 and the state ^3/2,1/2 of the n = 4 shell are almost degenerated for the chosen value 
of the deformation parameter d and consequently the occupation probability corresponding 
to one pair of protons is shared by the two states. 
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d d 



FIG. 2: Neutron (top) and proton (bottom) single-particle energy levels given in units of 
huj (= 41A- 1 / 3 MeV) and calculated with Eq.flZSD for 160 Yb (left) with X c = 74.9940 keV and 
162 Hf (right) with Xc = 78.2942 keV. The vertical lines indicate the single-particle configurations 
corresponding to the fitted deformation parameter d. The Fermi energy level resulting form the 
BCS calculations is also pointed out. 

The deformation parameter d affects both the single-particle and the collective degrees 
of freedom. Indeed, on one hand it is embedded in the strength dXc of the deformed mean 
field, and on the other hand it defines the energy of the core. One may therefore assert that 



16 




FIG. 3: The deformation parameter d (square) is presented as function of the nuclear deformation 
/?2. The fitted values of d can be interpolated by the straight line given by d = 19.55/32 — 1.49. 

the particle-core interaction induces a deformation effect on both the single-particle and the 
core motion. However, it can be easily checked that the ground band energies are not very 
sensible to the single-particle degrees of freedom. Indeed, the overwhelming contribution to 
the total energy of the ground band is due to the core because all the intruder particles are 
paired and do not carry any angular momentum. This fact implies that up to the first band 
crossing the whole angular momentum dependency is given by the core. 

Besides the deformation parameter d, the core energy is also parametrized by u)q and uj\, 
the strengths of the two boson terms. The core parameters d, Uq and oj\ are determined in 
the first approximation such that to reproduce the first yrast energy levels which are purely 
collective. The final value of the deformation d is fixed by achieving a consensus between 
the reproduction of the single-particle levels configuration and the best description of the 
angular momentum dependency of the total energy of the g-band up to the first backbending. 

The final touch to the formalism is made by fixing the strength C of the spin-spin interac- 
tion. The effect of the spin-spin interaction was presented in detail in Ref.j^ 1 . Basically, it 

[1] In Table II of Ref.f^j the values of C ■ 10 were listed. By a lamentable error the factor 10 accompanying 
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simulates the Coriolis force in the intrinsic reference frame and is actually the model Hamil- 
tonian term which is responsible for the pair breaking. Indeed, recalling the fact that the 
pair breaking is equivalent to the time-reversal symmetry breaking of the system it is then 
clear that it cannot be achieved by the qQ interaction and therefore the spin-spin interaction 
is necessarily demanded. It is found that this term does not have any effect on the energies 
of crossing bands up to the first critical angular momentum, but on the contrary has a strong 
effect on the moderate and high spin states in the yrast band. Because of this feature the 
strength C is fixed such that to reproduce the moderate and high spin yrast state energies. 

Apparently the number of parameters used is large, but three of them, d, G p and G n , 
are not freely changed when we pass from one nucleus to another. Indeed, the deformation 
parameter depends linearly on the nuclear deformation and therefore fixing it for one nucleus, 
for example by fitting the B(E2; + — > 2 + ) value, it is known for all remaining ones. This is 
shown in Fig. 3, where the fitted values of d were interpolated by a straight line. 

Also, the results for the strengths of proton and neutron pairing interactions can be 
interpolated by a function proportional with 1/A : 

Gv = MeV, B. = ^ MeV. (4.1) 

The A dependence of G T , t = p, n is quite close to the A-parametrization of the interaction 
strengths which interpolates the values obtained by fitting the even-odd mass difference: 

Gp = MeV, G. = ^ MeV. (4.2) 

B. Energies 

The energies of the rotational bands implied in the present model are approximated by 
the diagonal matrix elements of the model Hamiltonian between the projected states of 
the set (I2.12p and calculated using the parameters listed in Table I. The band mixing is 
achieved by orthogonalizing the projected states (I2.12p and then diagonalizing H (12. 2p in 
the resulting orthogonal basis. For a given total angular momentum J one obtains a set 
of four eigenvalues E™, with m = 1,2,3,4. The lowest energies E™ define the yrast band 
E{J). 



C was omitted. However, all results of the quoted reference correspond to the true values of C . 
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The band mixing is schematically shown in FigJH where all involved rotational bands and 
the resulting yrast band are plotted versus total angular momentum J. Similar dependence 
of the rotational bands on the angular momentum was obtained in Ref. 0, [l^] where the 
energies were computed only in a projected quasiparticle space with a relatively large number 
of single-particle states. As can be seen from Fig 0] the proton S'-band does not interact with 
the other bands or is very weakly interacting with the g-band at high spin states in the case 
of 162 Hf. Moreover, its energy is higher than that of other bands, such that it has no influence 
on the yrast band. Thus, the inclusion of this band is made for the sake of completeness, 
otherwise it could be ignored. However, the unperturbed proton S-band provides valuable 
information regarding the dynamics of the system's angular momenta. Indeed, the minimum 
displayed by both the neutron and proton S-bands in FigfJ] indicates the amount of angular 
momentum carried by the corresponding broken pair. This is suggested by the following 
reasoning. First of all one must note that the slopes of the curves from FigJD determine 
the rotational frequencies of the bands. The negative slopes of the neutron and proton 
S-bands at low spins imply a negative rotational frequency which is due to the core that 
must compensate the already high angular momentum realized by the decoupled broken 
pair. In the minimum point, where the slope vanishes, the core is no longer rotating and the 
total angular momentum is coming from the broken pair alone. Thus, the spin at which the 
S-bands show a minimum represents the angular momentum carried by the broken pair. 

Inspecting Fig]4]one finds that for all considered nuclei the neutron broken pair carries 
almost 8-10 units of angular momentum, while the angular momentum of the proton broken 
pair is about 6-8/?.. But as we already remarked, the second backbending is due to the 
crossing of the neutron S-band with the neutron-proton S-band and not with the proton 
one. Of course the 4qp band associated to two broken pairs, one of neutron and another of 
proton type has a different structure from a 2qp S'-band. As can be seen from FigJU such a 
band has an extended plateau which means that the total angular momentum is due to the 
both broken pairs without any core contribution. As a matter of fact the total spin where 
the plateau ends and the core starts to rotate is equal to the sum of the angular momenta 
provided by the broken pairs, which is around J = 16. 

For a better understanding of the multiple backbending phenomena, the theoretical re- 
sults and the experimental data are compared by means of backbending plots and the cor- 
responding energy spectra. The backbending plot is a graph which shows the dependence 
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FIG. 4: Energy trajectories implied in the band mixing are presented as function of the total 
angular momentum. The <?-band is represented by the straight line, neutron S-band by the dashed 
line and the proton S-band by the dotted line, while the dash-dotted line corresponds to the 
neutron-proton S-band. The yrast energies (circles) resulted from the diagonalization of the total 
Hamiltonian in the orthogonal basis (|2. 19|) are also visualized. 

of the moment of inertia on the angular frequency squared. If one adopts for the moment 
of inertia the following expression 

3= i^±6 (43) 

E(J+2)-E(jy K } 
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FIG. 5: Backbending plots for i&b Er, i58 Er, ibU Yb and iW Hf isotopes 



with experiment (circles). Experimental data are taken from Refs. 



comparing theory (squares) 



39rkj] 



where E(J) are the yrast energies, and defines the rotational frequency as 

MJ) = < ^dP~ W l [E{J + 2) ~ E{J)l (44) 
one readily obtains the experimental and theoretical backbending curves for the four nuclei 
treated here. These plots are shown in FigJS] where the description is limited to the experi- 
mental yrast states up to the spin J = 36 for 158 Er, 160 Yb and 162 Hf and J = 32 for 156 Er. 
The nature of states with angular momentum higher than 36 might be different from that 
of the states considered in the present work. Indeed, since the states density increases with 
the spin, one expects that a larger band admixture takes place. Even so, the number of 
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experimental states described here is enough to account for the most important features of 
the second moment of inertia anomaly. The smaller number of yrast states considered in 
the case of 156 Er is due to the fact that the states beyond J = 32 have not yet an angular 
momentum assigned. 

Coming back to the backbending plots of FiglSJ it is obvious that the double zigzag 
shape is reproduced quite well for all four nuclei. An especially good agreement is found for 
moderate spin states at the first backbending which is, indeed, very well reproduced in all 
cases. The second backbending is supposed to be less pronounced than the first one, because, 
as FigJH shows, the crossing angle between the neutron S'-band and the neutron-proton S- 
band is much smaller than the one between the g-band and the neutron S'-band. However, 
the experimental data offers a rather sharp second backbending for nuclei 156 Er, 160 Yb and 
162 Hf, while the theoretical calculations predict a smoother backbending behavior. In the 
case of 158 Er, the second observed moment of inertia anomaly is not a real backbending but 
a relatively weak up-bending. Note that, the theoretical results also predict an up-bending 
which is however much steeper. This is consistent with the results from FigJH where the 
crossing angle between the neutron and neutron-proton S-bands for this nucleus is very 
small. 

Concerning the comments on the pairing constants given in subsection B, the question 
which certainly arises is what is the effect on the backbending plot when we use the G T pa- 
rameters obtained by fitting the even-odd mass difference instead of the interaction strengths 
fixed as described above. The answer is given in FigEJ for 160 Yb. Indeed, comparing the re- 
sults corresponding to the two sets of G T , one may state that there is no significant difference 
between the two plots. Concluding, the real number of the free parameters is four. 

The good agreement between theoretical and experimental backbending plots is reflected 
also in the corresponding energy spectra. Thus, FigJT] suggests a very good agreement 
between the results of our calculations and the corresponding data, which is quantitatively 
expressed by relatively small root-mean-square (r.m.s.) values for deviations. Note that 
the energy spectra are better reproduced at high spins than at low spins, contrary to the 
backbending plots where the first backbending is better described than the second one. 
This happens because the backbending curves do not depend on the absolute energies of 
the angular momentum states, but on the energy difference between consecutive states 
and moreover through a quadratic law (hu) 2 which is more sensitive to small deviations. 
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FIG. 6: The backbending plot corresponding to two sets of parameters for the pairing strengths: 
one obtained by the fitting procedure described in the text (squares) and one fixed so that the 
even-odd mass difference (diamonds) be reproduced. 

Examining FigJTJ one remarks an increasing behavior of the critical energies with Z for the 
N = 90 isotones. This feature might be ascribed to the constant decrease of the deformation 
which increases the frequency of the collective rotation. 

The four nuclei treated here are 7-unstable. Thereby the collective motion of the N = 90 
isotones 160 Yb and 162 Hf can be well described by the 0(6) dynamic symmetry 431 ] . The 
softness of these nuclei points to a possible dynamic deformation which is increasing with the 
angular momentum. Indeed, judging by the behavior of the g-bands from FigJH the energy 
spectrum at lower spins is of the rotational type, while for larger spins it becomes more 
vibrational-like. This change in the energy spectra is most likely caused by the increase of 
the 7 deformation because the /3 is fixed for these nuclei. The structure of 156 Er is different. 
The observed collective spectrum of the 156 Er exhibits signatures of the E{5) dynamical 



symmetry [44| which is assigned to the critical point of the phase transition from the 0(6) 
to the U(5) symmetry. The critical point potential has a very extended minimum in the 
deformation parameter around the origin which corresponds to a spherical shape described 
by the U(5) dynamical symmetry. As a matter of fact, the observed nuclear deformation of 
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FIG. 7: Experimental and theoretical yrast spectra of Er, Yb, Hf and Er, with numerical 
values given in units of keV. The starting point of the backbendings are indicated for each nuclei by 
a dashed line. At the beginning of each spectrum one can find the corresponding root-mean-square 
(r.m.s.) values. 

156 Er is indeed small. In this case one can also have a variation of the (3 deformation along 
its flat minimum as the nucleus is increasing its rotation. 
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C. Angular momentum alignment 



In order to study the alignment of the angular momenta involved in the system's dynam- 
ics, it is useful to compute the averages of the involved angular momenta: 

J n (J n + l) = (<M|<«>, (4-5) 

J P (J P + 1) = (*j£\J*\*j£), (4.6) 

J f (J f + l) = «[|J7K[), (4.7) 

1(1 + 1) = (HZtflHZ)- (4-8) 

The deviation 

AJ= |j-(J c + J f )\, (4.9) 

is a measure for the departure from the full alignment of the fermionic and core angu- 
lar momenta, i.e. when J c + Jf equates the total angular momentum J of the system. 
All the average angular momenta (14. 5 1) - (14. 8 p and the deviation A J are plotted in Fig 18] as 
the functions of total angular momentum J. These plots reveal additional information for 
the backbending phenomenon. Indeed, from FigJH] one can extract the angular momentum 
carried by the neutron and proton broken pairs, the composition of the total angular mo- 
mentum, the critical spins of the band crossings, or one can even investigate the alignment 
of different angular momenta of the system. The difference between the values of the Jf, J n 
and J p before and after the critical angular momenta associated to the pair breaking, gives 
the amount of angular momentum carried by the broken pairs which is consistent to those 
determined from analyzing the plots of FigJH An interesting feature can be seen from FigJHl 
which is the essential difference between the two band crossings. Indeed, while the neutron 
angular momentum J n has a clear discontinuity reflected in a jump to a plateau of higher 
spin, the proton angular momentum has a steady increase extended around the critical an- 
gular momentum where the second band crossing actually takes place, although the curve 
changes substantially its slope. This was somehow expected due to the smaller crossing 
angle between the neutron and neutron-proton S'-bands. The smaller crossing angle means 
a larger range of the angular momentum where the bands are effectively interacting. The 
neutron and neutron-proton S'-bands start to interact from J = 22 for 156 Er and J = 24 for 
the rest of nuclei, and keep interacting afterwards. After this spin, the states are no longer of 
a pure nature and the nucleus is described by a coexistence of 2qp and 4qp states of broken 
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pairs. This is contrary to the case of the first band crossing where the interacting range is 
finite and very short, about 2 units of angular momentum. Investigating the behavior of the 
core anenilar momentum ,7„. one observes that it has a sudden fall at the first band crossing 
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FIG. 8: Color online. Expected values of the angular momenta corresponding to the neutron and 
proton broken pairs, the total fermionic angular momentum of the neutron and proton intruder 
orbitals as well as the core angular momentum. The deviation A J of the total angular momentum 
is also visualized. 



Concerning the angular momenta alignment, one remarks that before the band crossings 
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the alignment defect A J has a minimum and right after a local maximum. Note that here 
we deal with a rotational alignment and that is why the deviation A J decreases with total 
angular momentum. Even though, the full alignment A J = is not possible because of the 
fact that after the first band crossing the yrast states are of K ^ nature. However, at the 
beginning of the second band crossing, one finds that A J pa 0. This approximate alignment 
is due to the fact that the proton orbital starts to aid more consistently the fermionic angular 
momentum Jj when the neutron S'-band starts to interact with the neutron-proton S-band 
and the proton pair just slowly begins to break. This leads us to the conclusion that the 
angular momenta of the broken pairs first align to each other and only after that they align 
with the core angular momentum. The last alignment seems to be hindered, as shown in 
FigJH] where the angular momentum defect does not decrease after the second band crossing 
and moreover at some point it starts to increase in parallel with the core angular momentum. 
The increasing behavior of A J at high angular momentum states points to the fact that the 
rotation at high spins starts to work against the alignment between the core and fermionic 
angular momenta. 

D. Electric quadrupole transitions 

A very sensitive test of the wave functions describing the energy levels are the quadrupole 
transition probabilities. In Fig 19] one compares the numerical results provided by the formu- 
las from Sec. Ill with the corresponding experimental data available only for 156 Er, 158 Er and 
160 Yb. The parameters qi and (fe of the quadrupole transition operator are fixed by fitting 
the experimental B(E2) values and the obtained results are given in Table III. The theoret- 
ical and experimental values are also compared with the rotational limit of the quadrupole 
transition probability corresponding to the rigid rotor wave functions defined as: 

B(E2, J+ -> J' + ) rot = ^-Ql (c^)\ (4.10) 

with Qq fixed by fitting the first experimental transition probability B(E2,0 + — > 2 + ). The 
values of Qo corresponding to each considered nucleus are also given in the Table III. 

The transitions along the yrast band directly reflect the structural changes of the total 
wave function in the band crossing region. Indeed, investigating the theoretical points from 
Figj9] one notices that at the first band crossing only one transition is sizably hindered. 
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TABLE III: The results of the fitting procedure performed for the quadrupole transition proba- 
bilities shown in Figj9] are listed for each treated nucleus together with the Q$ value defining the 
values of B(E2) rot . 



Nucleus 


q x [W.u.] 1 / 2 


q 2 [W.u.]V2 


r.m.s. [W.u.] 


Qo [W.u.]V2 


i56 Er 


12.31060 


11.73900 


60.3801 


57.4669 


158 Er 


2.10613 


-1.58305 


74.4017 


77.34075 


160 Yb 


7.89159 


7.55747 


55.0754 


68.37170 



This indicates the fact that the interaction of the g-band with the neutron S"-band is weak 
such that the transition from Oqp to the 2qp nature is very sudden, taking place in the 
interval of no more than 2 units of total angular momentum. This behavior is also found 
in the experimental data, although in case of the 156 Er nucleus the minimum calculated 
transition is somehow shifted to the next transition in respect to experimental results. The 
situation at the second band crossing is essentially different because in this case both model 
states are of the quasiparticle nature which enforces the interband interaction leading to a 
less visible decrease of the transition probability with an extended minimum in the band 
crossing region. Looking at the experimental values, especially those before the first band 
crossing, we observe some significant deviations from the rigid rotor behavior. The largest 
deviations are obtained in case of the 156 Er nucleus. Judging by the moderately small values 
of the nuclear deformation @2 and of the obtained values for the deformation parameter 
d, it is not surprising that 156 Er deviates the most from the perfect rigid rotor case. The 
large discrepancy at the low spins between the experimental data and the predicted rigid 
rotor behavior could also be due to the fact that these nuclei are relatively sensitive to 
the shape fluctuations. This is, in fact, consistent with the previous comment about the 7 
softness of these nuclei. The oscillation of the transition probabilities before the first band 
crossing, although not yet well understood from the phenomenological point of view, it is 
well reproduced by the theoretical results. Indeed, even the unusual parabolic dependency 
on the angular momentum of the B(E2) values before the first band crossing in the 158 Er 
and 160 Yb nuclei is simulated quite well by the model predictions. 

Another feature which deserves attention consists of that the parameters q\ and qi for 
156 Er and 158 Er are quite different. One reason was already mentioned, namely that the 
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FIG. 9: Theoretical predictions for the reduced El transition probabilities are compared with 
experimentally available data for 156 Er, 158 Er and 160 Yb taken from Refs. 39h41||. The open symbols 



indicate experimental data with assumed or derived assignment and were not taken into account 
for the fitting procedure only in case of 156 Er nucleus. The rigid rotor limit of the B(E2) is also 
shown for comparison. 
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two isotopes have different deformation which makes 158 Er be closer to the rotor behavior. 
Another reason might be the fact that 156 Er reaches the conditions of a critical point of 
shape phase transition U(5) — > 0(6) exhibiting a E(5) symmetry which results in having 
a discontinuity in the strength parameters of the transition operator. Indeed, the ratios 
of the excitation energies E i +/E 2 + for the two isotopes are 2.315 and 2.743, respectively 
which have to be compared with the E(5) limit which amounts of 2.2. Therefore a smooth 
behavior of the q± and q 2 parameters in the isotopic chain of Er isotopes is expected to be 
broken at 158 Er which is close to the critical point of the shape phase transition. 



E. Gyromagnetic factor 



The magnetic dipole moment of the particle-core system is defined as: 



# = g c J c + 9fJf = gjJ, 



(4.11) 



where g c and g$ denote the gyromagnetic factors of the core and fermionic subsystems, 
respectively. The structure of the total wave function is reflected by the total gyromagnetic 
factor gy. 



9j = 9c + 



9f~9c 



J f (J f + l)-J c (J c + l) 
J(J+1) 



For the core gyromagnetic factor one takes the rotational value 



(4.12) 



(4.13) 



given in units of nuclear magneton /xjv, where Z c and A c are the nuclear charge and the 
mass number of the core: 



T irh-i 



pair 



(4.14) 
(4.15) 



with the expected number of neutron and proton pairs determined from the BCS equations 
and given in Table II. 

As for the fermionic gyromagnetic factor, it is obtained from the following decomposition 
of the fermionic magnetic moment: 



fif ~ 9fJf — 9nJn + 9pJp, 



(4.16) 
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FIG. 10: Color online. Calculated gyromagnetic factor for yrast states (circles) given in units of 
nuclear magneton is represented as function of angular momentum. There are also visualized few 



experimental values (squares) taken from Refs. 
the gyromagnetic factor. 



391442] together with the rotational limit Z/A of 



which gives an expression for g$ in terms of J n , J p and Jf similar to (I4.12p . 



9f = 9 P + 



9n ~ 9p 



1 + 



J n (J n + l) - J p {J p + l) 



(4.17) 



Jf(Jf + 1) 

Knowing that the intruder neutrons are from the iiyz orbital, and the intruder protons 
are from the hu/2 orbital, one obtains the following values for the proton and neutron 
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gyromagnetic factors 



9n = gl n) + (gi n) - gl n) )/is = f| = -0.22^, (4.18) 

9v = 9? + (# - rf P) )/H = 1-29^. (4.19) 

ree values were 



For the above calculation we used the free value of the gi while for g s the 
quenched by the factor 0.75, which accounts for the nuclear medium effect 



45|, 



g? = 0, gf = lfx N , ^ = -3.8256 x 0.75^, ^ = 5.5855 x 0.75/^. (4.20) 

The total gyromagnetic factor is plotted in FigJlO] as function of the total angular mo- 
mentum J. Its change in the behavior reflects the transition from states of different nature. 
Before the first band crossing its value is almost constant and close to the rotational limit, 
although slightly overestimated. Of course, even if the nature of the g-band is collective, it 
is far from being perfectly rotational as it is suggested by the small values of the deforma- 
tions d and (3 2 from Table I. Indeed, it can be seen from Fig{10] that the departure of the 
gyromagnetic factor from its rotational limit Z/A before the first band crossing is bigger 
for 156 Er and 162 Hf nuclei, which turn out to be the less deformed ones. At the first band 
crossing the gyromagnetic factor has a sudden fall down, reaching very small values where 
the total magnetic moment almost vanishes. This discontinuity marks the change of the 
yrast band from Oqp to a 2qp neutron character. The fall of gj at the first band crossing 
is due to the negative value of the neutron gyromagnetic factor coming from the decoupled 
neutron pair. After the first band crossing the rotation of the core starts to dominate and 
the gyromagnetic factor increases almost linearly with J. This trend keeps up to the second 
band crossing where the ascendant slope becomes bigger due to the positive value of the 
proton gyromagnetic factor coming from the proton broken pair. The second band crossing 
is reflected in an inflexion point of gj as function of J. This is consistent to the slowness 
of the consequent breaking of the proton pair which does not offer a jump like in the case 
of the first band crossing. The growth of the gj persists only for a few states and then it 
comes to a saturation in the vicinity of the rotational limit value. As a matter of fact, the 
mentioned plateau begins at the spin where the second backbending ends. Few remarks are 
necessary regarding the comparison of calculation results with the experimental values of the 
gyromagnetic factor. Leaving aside the nuclei 156 Er and 162 Hf where relevant experimental 
data are lacking, the other two reproduce quite well the sudden fall of gj at the first band 
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crossing. An especially good agreement between theory and experiment is obtained for 158 Er 
where not only the discontinuity of the gyromagnetic factor but also its absolute values are 
reproduced. 

Before closing this section, we would like to comment on the obtained values of some 
of the model parameters. First of all, one notes the linear dependence of the deformation 
parameter d on the nuclear deformation fa. This property can be used to approximately 
determine the deformation d for other nuclei from the rare earth region. The numerical 
values of the deformation parameter d are in the range of values determined in Refs. 46, 4?| 
for other isotopes of the nuclei treated in this paper. This feature pleads in favor of both the 
CSM formalism and the present approach. The other parameter which deserves a special 
attention is the strength of the spin-spin interaction. Although such an interaction was 
already used in connection to the backbending phenomena 2l|, here it brings an essentially 



different contribution. First of all in Ref. 



21], the spin-spin interaction was found to be 



repulsive while in the present model it can be both attractive and repulsive. Indeed, the 
spin-spin interaction matrix elements are going from negative to positive values in the 2qp 
and 4qp bands as well as in the corresponding non-diagonal matrix elements. The picture 
is opposite for negative values of the strength C, as happens in the case of 162 Hf. It is 
interesting to mention that the second backbending in 162 Hf, is difficult to explain due to its 
unexpected sharpness. Indeed, before the second backbending was experimentally observed 
in 162 Hf, the CHFB calculations predicted for this nucleus a small up-bending or even no 

n 

backbending [48|. As a matter of fact in our approach the reproduction of the second 
backbending in this nucleus was possible only by choosing a negative value for the spin-spin 
strength C . This feature proves the importance of the spin-spin interaction in explaining 
the backbending phenomenon which thus appears to be the result of an interplay between 
the Coriolis-like force and the Qq interaction. 



V. CONCLUSIONS 



The present model provides a consistent explanation for the pair breaking process in 
connection to the rotational alignment of the angular momenta involved in the system. 
Using simple arguments one determines the critical angular momentum J where the pair 
breaking takes place. While the neutron pair breaking takes place at J = 10 or 12, one 
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cannot accurately say at what angular momentum the proton pair is broken because at high 
spin states the crossing bands interact within a larger interval. This is suggesting that the 
proton pair breaking is a slower process than the neutron pair breaking. 

Concerning the rotational alignment, it is found that the proton and neutron angular 
momenta first align to each other and only after that they align to the core angular mo- 
mentum. The full alignment between the fermionic and the core angular momenta cannot 
be achieved due to the intrinsic properties of the higher spin states which are of the K — 1 
and K = 2 nature. However, strong alignments are obtained at the band crossing critical 
angular momenta. Another interesting result of the present approach is that the rotational 
alignment lessens after the second backbending, which is pointing to the fact that the 2qp 
and 4qp bands still interact even after the band crossing. 

The first backbending manifests itself in the gyromagnetic factor plot by a big fall down 
of gj. By contradistinction, the second backbending is reflected by an inflexion point in the 
above mentioned plot. 

The effect brought by each term of the model Hamiltonian on the spectrum in the region 
of the band crossing is in extenso analyzed. In this way the free parameters acquire a well 
established significance. 

Along the time, various versions of angular momentum projection has been used with the 



aim of describing the backbending phenomena |8Hl2 



14 



18 



421]. 



What distinguishes our model from the others? First of all the three components of 
neutrons, protons and the core are described by deformed wave functions. Moreover, the 
mean fields of neutrons as well of protons are derived from the particle-core coupling term and 
thereby the three components have similar deformation properties. The total wave function 
describing the nucleus in the laboratory reference frame is obtained by angular momentum 
projection procedure from the product of the mentioned three deformed functions which, is 
not an easy task. We suspect that due to the specific construction, the wave function has 
a complex structure which allows to describe quantitatively the spectra in the region of the 
two backbendings. The accuracy of description is reflected not only in the backbending plot 
but also by transition probabilities (Fig.9) and gyromagnetic factors (Fig. 10). 

Note that the core is described by projecting out the angular momentum from a coherent 
state and by an anharmonic boson Hamiltonian. Therefore the core moment of inertia is not 
constant but depending on the angular momentum. In this context we could assert that our 
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model is on a par with those particle-core approaches using a variable moment of inertia. 

As a final conclusion one can say that the present formalism is able to describe quantita- 
tively the double backbending phenomenon. Moreover, a consistent qualitative explanation 
of the combined contribution of the pair breaking and rotational alignment to the backbend- 
ing phenomenon is provided. 
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